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Abstract 

We use the computer algebra system GRTensorll to examine invari- 
ants polynomial in the Riemann tensor for class B warped product space- 
times - those which can be decomposed into the coupled product of two 
2-dimensional spaces, one Lorentzian and one Riemannian, subject to the 
separability of the coupling: 

ds 2 = dsl t (u, v) + C(s 7 ) 2 ds| 2 (6, cf>) (1) 

with C{x'<) 2 = r(u,v) 2 w(8,<f)) 2 and siff(Ei) = 0, sig(E 2 ) = 2e ( e = 
±1) for class Bi spacetimes and sig(E\) = 2e, sig(Y>2) = for class 
£>2. Although very special, these spaces include many of interest, for 
example, all spherical, plane, and hyperbolic spacetimes. The first two 
Ricci invariants along with the Ricci scalar and the real component of 
the second Weyl invariant J alone are shown to constitute the largest 
independent set of invariants to degree five for this class. Explicit syzygies 
are given for other invariants up to this degree. It is argued that this set 
constitutes the largest functionally independent set to any degree for this 
class, and some physical consequences of the syzygies are explored. 



1 Introduction 

It is natural to suspect that invariants constructed from the Riemann tensor via 
contractions^ should play some important role in general relativity. However, 
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it is only fair to say that this importance has yet to be fully realized. For 
example, it has been known for a very long time that this type of invariant fails 
to distinguish inequivalent spacetimes . Indeed, the inclusion of derivatives of 
the Riemann tensor does not overcome this shortcoming || . The inequivalence 
problem requires a more sophisticated approach (e.g. pj). Invariants are of some 
use in the study of spacetime singularities (distinguishing the scalar polynomial 
(sp) type (e.g. [p|) which is equivalent to a curvature pp type Q), but even 
here they fail to distinguish shell crossings from more significant singularities. 

Despite these shortcomings, the problem of finding a complete set of invari- 
ants for spacetime has attracted considerable recent effort [0, ||, ||, |l(J . Here we 
set out to study spacetimes which are the product of two 2-dimensional spaces, 
one Lorentzian and one Ricmannian, subject to a separability condition on the 
function which couples the 2-spaces. That is, we consider metrics of the form 

ds 2 = ds^ (u, v) + C(x^) 2 dsl 2 (0, 0) (2) 

subject to the restriction 

C(x"<) 2 = r(u,vfw(e,4>) 2 . (3) 

The spaces are known as warped products of class B [[fl], with sig(Y,\) = 0, 
sig(T,2) — 2e ( e = ±1) for class B\ and sig(Ei) — 2e, sig(£ 2 ) = for class B 2 . 
A metric of sufficient generality for class Bi is given by^ (e.g. |Q]) 

ds 2 = -2f(u, v)dudv + r(u, v) 2 g(9, (j)) 2 (de 2 + d<p 2 ) (4) 

and for class B 2 by 

ds 2 = f(u, v) 2 (du 2 + dv 2 ) - 2r(u, v) 2 g(9, <j))d6d<j). (5) 

There are many examples of class Bi spacetimes (for example, all spherical, 
plane, and hyperbolic spacetimes). Examples of spacetimes of class B2 are 
considerably harder to find. It is known that the only physically interesting 
energy-momentum types of class B2 spacetimes are non-null electromagnetic, 
A-term, or vacuum |13J . 



2 The number of invariants required 

For a general Petrov type D metric, we can align the tetrad along the principle 
null directions of the Weyl tensor to find a frame in which there are in general 
12 non-zero components of the curvature (complex ^2, the Ricci spinor com- 
ponents, and the Ricci scalar). Since this alignment leaves a two parameter 
group of rotational freedom in the frame ('spin' and 'boost'), the number of 
independent invariants for a type D spacetime is reduced by this dimension to 
10. 

2 We assume that all functions are twice continuously differentiable. We do not consider 
surface layers. Note that we have chosen to write S2 in conformally flat form for class B\ . 
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Specializing to the given class of spacetimes, we find that if the above frame 
alignment is carried out, then a number of the Ricci spinor components are also 
reduced to zero. For class B\ we are left with the real components $00, $11, 
and $22, while for class -B2, only <E>n and $02 = $20 remain non-zero. The 
number of independent degrees of freedom arising from the Ricci components 
is thus reduced to 3. However, one of these degrees can be removed by making 
use of either the spin (for B\) or boost (for B2) freedoms still remaining in the 
frame. 

The number of independent invariants that we should expect to find for the 
given class of spacetime is thus 

2 + 2 + 1-1 = 4, 

corresponding respectively to the Weyl and Ricci spinor freedoms, the Ricci 
scalar, and the dimension of the invariance group of these spinors. 

An example of how the details of this calculation can be carried out using 
GRTensorll is given in Appendix A. 



3 The invariants 

The invariants we begin with here are the sixteen suggested by Carminati and 
McLenaghan augmented with an invariant equivalent to M suggested recently 
by Zakhary and Mcintosh ||.^] We adopt the definitions in [pj for the setQ 

CM = {R,r 1 ,r 2 ,r 3 ,wi,W2,mi,m2,m3,m4,m 5 } (7) 

and augment these with 

ZM = m 6 = ^ ABC d^ AeAb ^ B ecd ^ C fAb ^ DFCD . (8) 

In total this set constitutes six real and six complex invariants. The complete- 
ness of this set has been examined recently || . It is certainly complete to degree 
five (degree being the total number of Weyl, conjugate Weyl, and Ricci spinors 

3 As an aside we note here that if we consider, in the usual way (e.g. Jl^t), the secular 
equation \S a t — Xg a t\ = 0, where S a t is the trace-free Ricci tensor, and g a i, is the metric 
tensor, the characteristic polynomial which follows is ctiX 4 + oj2^ 3 + 03A 2 + ct^X + 05 = 0. 
Since the Oj s are derived from the components of S a b and g a f, by algebraic operations, they 
can be called 'algebraic' invariants of the tensor S a b- It can be shown that for this class of 
spacetimes the following relationships hold between the invariants ri, rg, rz and ai, c?3, 04, 05 
(02 = is a consequence of S a b being trace-free): 

n=--23- r9 = Mi. r _ _2l _ *s_ (a) 

11 2oti' 12 Bar' ' 3 JS( 4a 1 ' [ -> 

4 Two scalars found frequently in the physical literature are the Kretschmann scalar 
(RiemSq = R a bcdR abcd \ Rabcd the Riemann tensor) and the square of the Ricci tensor 
(RicciSq = R a tR ab , R a b the Ricci tensor). In terms of the set (pf) these are given simply 
by RiemSq = R 2 /6 + 8(n + Re(w\)) and RicciSq = (R 2 + n)/4. 
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whose contractions make up an invariant )[]. In what follows we explicitly reduce 
this set of eighteen invariants to four for the spacetimes under consideration by 
way of syzygies jl5| . All syzygies have been developed and tested with the aid 
of the computer algebra package GRTensorll |l6|].p| 



4 Syzygies 

Syzygies follow from a knowledge of which spinor components are non-zero. In 
the present case the important features are that these components are real and 
that the Weyl and Plebanski tensors (e.g. Jig] ) share the same principal null 
directions and are of type D (or O). Since spacetimes of the form (2) with (3) 
are of Petrov type D (or O), we have the well-known elementary syzygy^ 

6wj -wl = 0. (10) 

Furthermore, for any Petrov type D (or O) spacetime we have |l7j 

(3m2 — w\ri)wi — 3toiW2 = (11) 



and 

(3rri5 — w\mx)w\ — 3m 3 w2 — 0, (12) 



5 It is a simple (though tedious, unless carried out electronically) matter to write down 
expressions for all index contractions among any set of n curvature spinors. If this is done to 
degree five, it is found that for general spacetimes every invariant which is not included in the 
CM+ZM set satisfies an identity allowing it to be written as a rational function of members 
of the set. Thus, there are no invariants of degree less than or equal to five to be added to 
the CM+ZM set. [|u| 

B GRTensorll is a package which runs within MapleV. It is entirely distinct 
from packages distributed with MapleV and must be obtained independently. The 



(j-Klensorli software and dr 
Web from the address 


cumnntatioTi is rlistri huted trenlv nn the World-\ 
tittD: / /www. astro . queensu. ca/~g"rt elisor /GRHome .html 


ittp : //www. maths . soton. ac . uk/~dp/grtensor/|. Worksheets which reproduce the 



or 



7 A11 syzygies can in fact be obtained algorithmically. By appropriately defining the "de- 
gree" of a syzygy, the most general syzygy can be constructed for any particular "degree". 
Expanding and collecting like terms allows the coefficients to form a linear system of equations 
which are then solved to find the coefficients of the syzygy. If there are no nonzero solutions af- 
ter the linear system is solved, then the syzygy must be of a different "degree" . The definition 
of "degree" in this case imposes an ordering, allowing a sequential search to be executed. In 
any particular case (e.g. as considered in this paper) the number of degrees of freedom can be 
counted and so the search for syzygies can be made exhaustive. Whereas degree is commonly 
used as a single parameter, as given in Section 3 above, a multi-parameter definition of degree 
is more convenient here. For example, an invariant can be given a "degree" [A,B,C], where 
A refers to the exponent on the Weyl spinor, B refers to the number of multiplications in the 
Ricci spinor, and C refers to the exponent on the conjugate Weyl spinor. In this notation, for 
the types of spacetimes we consider here, only two parameters are necessary, [A,f?]. Given a 
rule of composition: for invariants II = [Al, Bl] and 12 = [A2, B2], the degree of II * 12 is 

[Al, Bl] + [A2, B2] = [Al + A2, Bl + B2], (9) 

the parameter space given by [A, B] can then be searched exhaustively for the existence of 
syzygies. 
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where mi is the complex conjugate of mi.Q Equivalent syzygies have been given 
by Zakhary and Mcintosh for their set of invariants P] . We observe the following 
additional syzygies p]| 

6m,4 + wir2 = (13) 

and 

m 3 - m 2 = 0. (14) 

Further, we observe that Jl^] f\ 

(-12r 3 + 7rl)w imi - (12r| - 36rir 3 + 17r?)w 2 = 0, (15) 

2(3m-6 — miri)w2 + m^wi — 0, (16) 

and that 

(-12r 3 + 7r\f - (\2r\ - 36rir 3 + Ylr\f = 0. (17) 

The syzygy ( p7| ) is the counterpart for the Plebanski tensor to the syzygy 
(0) for the Weyl tensor]^ We expect that the syzygies (|l0|) through (|l7|) will 
hold in a wider class of spacetimes than type B warped products. 



5 The independent set of invariants 

The syzygies ([ll]) through ( |l7| ) remove the independence of m,2, m-5, m^, m 3 , 
mi, rriQ, and r 3 respectively. According to the syzygy ( |Io| ) the sign of W2 (which 
changes with signature) cannot be obtained from w\ whereas w\ > 0. This 
suggests that for metrics of the form (2) with (3) we take the set 

{R, ri ,r 2 ,W2} (18) 

as the independent set of scalar polynomial invariants, satisfying the number of 
degrees of freedom in the curvature. 

We note, however, that these do not necessarily constitute a complete set in 
the sense of classical invariant theory, where the set is considered complete only 
if every invariant not contained in the set can be written as an integral rational 
function of the other members of the set ^TJ| . This is clearly not the case here, 
for wi = ^ can n °t be written as an integral rational function of W2 = ^l: nor 
vise versa, whereas r 3 is determined as the root of a cubic from 

8 For type D w\ and W2 are non-zero, whereas for type O (or TV or III) both w\ and W2 
vanish. In the set CM, for a general Ricci tensor, there remain 10 invariants for type D, 9 for 
type III (since miR=m2C=mc l R=0 for the real (R) and imaginary (C) parts), 7 for type N 
(since mi C=m,o=m,^=p) and of course 4 for type O. 

^Regarding syzygy ( |15[ ) we note that mi = => 7rJ = 12r^ and that 7rf = 12r$ mi = 
or miuii = 3ti)2ri. 

10 We thank a referee for pointing this out. 
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From a practical standpoint, however, we are more concerned with isolating 
the functional degrees of freedom in the curvature, since by considering a more 
general interdependence (ie. by allowing non-integer exponents in our polynomi- 
als) we can in principle solve for any other of the remaining invariants in terms 
of the given set. By inspection of the invariants listed in (|l^), it is clear that 
they are functionally independent and that they make up the required number 
of degrees of freedom (i.e. 4) for the class of spacetimes in question, and in this 
sense we can consider them to be a 'complete' list. 

The invariants are reproduced in Appendix B. 



6 Some consequences of the syzygies 

The syzygies are certainly of value for computational efficiency. However, under- 
standing the physical benefit of the syzygies is somewhat more difficult. Mixed 
invariants have been used to provide a number of alignment theorems for the 
perfect fluid and Maxwell cases j7| and to classify the Riemann spinor |H| . Our 
purpose here is to provide some physical consequences of the syzygies presented 
by way of expressing them in terms of the electric and magnetic components 
of the Weyl tensor, and to show how the syzygy ( |l7| ) can be used to impose 
restrictions on physical quantities in the energy-momentum tensor. 

It is convenient to revert to classical notation. First we consider the Weyl 
tensor (C a bcd)- The "electric" and "magnetic" parts are defined, as usual, by 

E ac = C abcd u b u d , (19) 

and 

Hac = C* abcd u b u d , (20) 

where u a is a unit timelike vector and C* bcd = \r\abejC^ c d the dual tensor. 
Writing 

Qab = E ab + iH ab (21) 

and 

C abed — C a bcd + *C a bcd) (22) 

directly from the expression 

4v,[ a Q b ][ d U c ] + 9a[cQd]b ~ 9b[cQd]a + ir]abefU e U[ c Q d ] f + ir] c defU e U[ a Qb] f (23) 

for -\Cabcd (e.g., 10) it follows that 

CabcdC — QabQ ab (24) 



and that 



C ab cd C e fC c d } = Q a b Q b c Q c a . (25) 
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Equivalently, it follows that w\ (= 21) is given by 



^-(E ab E ab - H ab H ab ) + -(H ab E ab ) (26) 
lb o 

and that W2 (= 6 J) is given by 

1 7 

/o rpa rrb ttc rpa rpb rpc \ , / rja rrb ttc q rpa rpb ttc \ /97\ 

— (dti h tt c tl a — & b H, c £j a) + 7^{-H b n c-H a — (,£) c -tt a ). (Z ( ) 



Relations (26) and ( pTj ) hold m general and are independent of the timelike 
vector chosen to split the Weyl tensor into its electric and magnetic parts. 
These expressions can, for example, be used to clarify the relationship between 
the criterion of Bel and those of Misra and Singh. The Bel condition |2^] is that 
W\ = W2 = if and only if the spacetime is of Petrov type N or III (or, of course 
O). The condition of Misra and Singh || can be stated as E ab E ab = H ab H ab 
and E a b E b c E c a = H a b H b c H c a . Expressions (§|) and @ show that these 
conditions are inequivalent, contrary to the claim by Zakharov p4| ] . 

Independence from the timelike vector chosen to split the Weyl tensor does 
not, of course, extend into the mixed invariants. For example, it follows that 
the real and imaginary parts of 8mi are given by 

C abcd S ac S bd = 2EZS b S c a - 2E ab S ab S cd u c u d + S a c u c E ab S b d u d , (28) 

and 

C* abcd S ac S bd = 2H£S b S c a - 2H ab S ab S cd u c u d + S a c u c H ab S b d u d , (29) 

respectively, where S ab is the trace-free Ricci tensor. Expressions for the next 
degree of mixed invariants (rri2 in the notation of ) have been evaluated but the 
results are rather unwieldy (they contain 16 independent terms). The important 
point is that with the appearance of the vector S^U° and scalar S cd u c u d it 
becomes necessary to introduce an explicit form of S ab in order to proceed. 
This avenue is explored briefly below. 

Relations equivalent to (^4j) (or (|26| )) have been reexamined recently by 
Bonno r pE| . Relations equivalent to (|25|) (or (|27|)) have been given by Mcintosh 
et al. |26| for Petrov types I (or D). 

Specializing now to real \E , 2 we observe that 

H ab E ab = 0, (30) 

and that 

Tja rrb ttc o rpa Tpb TTC /qi\ 

tl b tl c tl a — 6tj b £j c tl a . 

Moreover, from (9) we obtain the further relation 

24(3E\H b c H c a - E a b E b c E c a f = {E ab E ab - H ab H ab ) 3 . (32) 
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It is perhaps worth emphasizing the fact that relations ( p0| ) , ( |3l| ) , and ( |32| ) hold 
for all Petrov type D spacetimes with real \&2< This includes a much wider class 
than class B warped product spaces. 

We now wish to explore any restrictions that the syzygy ( p"7| ) might impose on 
physical quantities in the energy-momentum tensor. First, however, recall that 
at most three Ricci invariants are independent. Finding syzygies associated with 
Ricci invariants of higher degree^] is a purely algorithmic procedure.^ Ricci 
syzygies up to rio are given in Appendix C. 

To begin, consider a perfect fluid, 

Tab = (fi + p)u a u b + pg a b (33) 

where p is the energy density, p is the (isotropic) pressure, and u a is the nor- 
malized timelike flow vector. The Ricci invariants in this case all reduce to the 
form 

r n = m n (p + p) n+1 , (34) 



where m n is a constant. The syzygy (17) reduces to the identity = 0. In 



creasing the complexity of the energy-momentum tensor slightly does not offer 
anything new. For both the case of anisotropic pressure p^] 

Tab = (p + P2)u a u b + (pi - p 2 )n a n b + p 2 g a b, (35) 

where n a is a normalized spacelike vector orthogonal to u a ( u a n a = 0), and for 
the case of a perfect fluid with a heat flux |2i| , 

Tab = (p + p)u a u b + pg a b + 2u( a <? 6 ) , (36) 

where q a is the heat flux vector such that u a q a = 0, the syzygy (|l7]) again 
returns the identity = 0. No information on possible relationships between 
the physical quantities that describe the spacetime can be gained from the 
syzygy ( |l7j ) in these cases. To obtain any restriction, a more complex physical 
decomposition must be considered. The simplest energy-momentum tensor for 
which the syzygy provides any restriction is one with anisotropic pressure and 
a heat flux: 

Tab = {p + P2)u a u h + (pi -p2)n a nb + P29ab + 2u( a g fc) . (37) 
In this case the syzygy ( ]Tt| ) reduces to 

(q a q a - {n a q a ) 2 ) 2 {pl - p2) 2 P = 0, (38) 



11 These are constructed analogously to n, r2, and . For example, 

r 4 = -±S a b S b c S c d S d e S e a = ^b.ab^C^^S^^B^b*^ 

12 All possible combinations of ri, and and r n (for each n > 4) are combined at degree 
n + 1. Expressing this combination as a polynomial in the Ricci components, the coefficients 
of any one of the Ricci components will form a linear system of equations that can be solved 
to find the numerical coefficients of the syzygy. 
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where P is a polynomial of degree 4 in 



{p,Pi,P2,n a q a } (39) 

and degree 3 in q a q a . If pi ^ p2 and q a q a ^ (n a q a ) 2 then the syzygy ( p7[ ) 
reduces to the equation P = and this reduces the number of independent 
physical scalars to four, the number of degrees of freedom. P is given by 

-p2 2 p A - pi V + 2plp2p 4 - 2p2V - Ap2(n a q a ) 2 p 3 + A{n a q a ) 2 plp 3 + 
2q a q a p2p 3 - 2q a q a p\p 3 - 2plV + 2pl 2 p2p 3 + 2p\p2 2 p 3 - (q a q a ) 2 p 2 - 
2p\ 3 p2p 2 + 6(n a q a ) 2 pl 2 p 2 - pl 4 p 2 + 8q a q a p2 2 p 2 + 2q a q a p\ 2 p 2 - 2p\p2 3 p 2 - 
6p2 2 (n a q a ) 2 p 2 - p2 A p 2 - I0q a q a plp2p 2 + 6pl 2 p2 2 p 2 + 8(q a q a ) 2 plp - 
30p2 2 (n a q a ) 2 plp + 2q a q a p2 3 p + 6p2 3 (n a q a ) 2 p + 8q a q a pl 3 p + 2p\ 3 p2 2 p + 
30pl 2 (n a q a ) 2 p2p + 2p2 3 p\ 2 p + \8q a q a p2{n a q a ) 2 p - Qpl 3 {n a q a ) 2 p - 
2plp2 A p + 10q a q a p2 2 plp - W(q a q a ) 2 p2p ~ 2p2pl A p ~ I8q a q a (n a q a ) 2 plp - 
20pl 2 q a q a p2p + 2pl 2 p2 2 q a q a - I8q a q a p2 2 (n a q a ) 2 + Wpl 3 (n a q a ) 2 p2 - 
I0p2 3 (n a q a ) 2 pl - p2 i pl 2 - (q a q a ) 2 p2 2 - 4(n a q a ) 2 pl i - 8pl 3 q a q a p2 + 
4q a q a pl 4 - 8pl(q a q a ) 2 p2 + 4{q a q a ) 3 + 2p2 3 plq a q a + 2p2 3 p\ 3 + 
54(n a q a ) 2 plq a q a p2 + 4p2 i (n a q a ) 2 - M(n a q a ) A plp2 + 27p2 2 (n a q a ) i + 
27(n a q a ) i pl 2 - mq a q a (n a q a ) 2 p\ 2 - pl 4 p2 2 + 8(q a q a ) 2 pl 2 . 
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7 Appendix A: Curvature components of the 
Class B\ warped product spacetimes 

The following session demonstrates how GRTensorll can be used to hx the frame 
components of the curvature relative to a 'standard' form by applying SX(2, C) 
rotations of the frame vectors. The standard frame is defined so as to make use of 
this rotational freedom in order to reduce the number of independent functions 
among the curvature components. Once the frame is fixed as far as possible, 
the subgroup of SL{2, C) under which the resulting curvature components are 
invariant is reported. The curvature components in the standard frame are 
called 'Karlhede invariants' for their usefulness in determining the equivalence 
of spacetimes (e.g. 0). 

Within MapleV, we restart a fresh session and initialize the GRTensorll 
package. 

> restart: 

> grtwQ: 

> qload (Bl) ; 
Calculated ds for Bl 

Default spacetime = Bl 

Coordinates 
x a = [u, v, 6, <f>] 
Line element 

ds 2 = -2f(u,v) du d v + r(u, vf g(0, <j,f dd 2 + r(u, vf g(0, (f>) 2 d(f> 2 

We have loaded a line element corresponding to the class B\ warped prod- 
uct spacetimes. However, we will eventually wish to examine the spinor 
components of the curvature. The function nptetradO can be used to 
automatically create a null (specifically, Newman-Penrose) tetrad from the 
metric. This tetrad will then be used to calculate the curvature spinors in a 
spin frame corresponding to the null frame. 

> nptetradO ; 

The metric signature of the Bl spacetime is +2. 

In order to create an NP-tetrad, the signature of g(dn,dn) will be 
changed to -2. 

Continue? (l=yes [default] , other=no) : 

> l; 

The signature of the Bl spacetime is now -2. 
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For the Bl spacetime : 



Basis inner product 



1 











.(«)(») 



1 











r, 











1 







1 







Null tetrad (covariant components) 



l a = [0, f(u, v), 0, 0] 



n a = [1, 0, 0, 0] 



m a = 0,0 



', ~lV2g(6, 4>)r{u, v), ~ V2g(6, <jy)r(u, v) 




The null tetrad has been stored as e(bdn,dn). 

The curvature spinors are denned as the NP scalars {WeylSc, RicciSc, and 
Lambda} in the standard GRTensorll package. However, they are also de- 
fined specifically as spinors {WeylSp, RicciSp, and Lambda} in the spinor .m 
package, which we now load. By calculating the latter objects, we can make 
use of the spinor packages' facilities for performing frame rotations and fixing 
a standard form for symmetric spinors. 

> grlib (spinor) ; 

Spinors and spacetime classification routines 
Version 0.9 4 Dec 1997 

> greale (WeylSp , RicciSp , Lambda) ; 

Basis/tetrad related object definitions 
Last modified 5 February 1997 

Calculated WeylSp for Bl 

Calculated RicciSp for Bl 

Calculated Lambda for Bl 
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The next command expands the tensor components (the first argument, 
tells the command to refer to the previously calculated spinors) and applies a 
function which aliases the output to make it more readable. 

> gralter (_ , expand , autoAlias) ; 

Component simplification of a GRTensorll object: 

Applying routine expand to object WeylSp 

Applying routine expand to object RicciSp 

Applying routine expand to object Lambda 

Applying routine autoAlias to object WeylSp 

Applying routine autoAlias to object RicciSp 

Applying routine autoAlias to object Lambda 

> grdisplay(_) ; 

For the Bl spacetime : 

*2 = 

1 9<t> 2 1 99 2 1 1 fu, v 1 fufv 1 9e,e 1 9 

' <-i p ^ Or s. r O rQ ' f . O Q ~T ✓» ' 



6r 2 g 4 6r 2 g 4 3 r f 3 r 2f 6 p 1 6 / 3 6r 2 g 3 6 r 2 

Tu fu 7*u, u 



rf 



1 r u r v 1 g 2 1 ff 9 2 _ 1 1 /„ /„ 1 50, 6 1 

^1 1 O r * a O A ' a OA a rO ' a 



2 r 2 J 4 r 2 g 4 4 r 2 ff 4 4 /2 1 4 J3 4 r 2 ff 3 4 r 2 ff 3 



2 2 



r / 2 r / 3 



1 -1 ^ OA*** o r * ^ /*~r' 



12 r 2 g 4 12 r 2 g 4 6 r 2 f 3 rf 12 / 2 12 / 3 12 r 2 3 3 

_j_ g^,0 
12 r 2 g 3 

At this point, we see that the Weyl spinor is already in what we would regard 
as a 'standard' form. That is, the I and n vectors of the frame are aligned along 
the principle null directions so that the only remaining Weyl component is ^2- 
In this form, the Weyl spinor is invariant under the SL(2, C) subgroups of spins 
and boosts, as we can check using the function isotestO. Thus by requiring 
the Weyl spinor to take this form, we have already lost four degrees of frame 
rotational freedom. 

> isotest (WeylSp) ; 

{Inswap, Spin, Boost} 
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(inswap refers to the discrete group corresponding to a swap of the I and n 
basis vectors). We can apply the same procedure to the Ricci spinor, 

> isotest (RicciSp) ; 

{Inswap, Spin} 

and learn that the Ricci spinor is, in fact, not invariant under boosts of 
the frame. That is, we can remove this freedom from the frame by using 
it to restrict the components of the Ricci tensor which are non-invariant 
under boosts so that they fit some standard form. A set of standard forms 
for symmetric spinors is known to the function dytrgenO, which can now 
be used to generate the appropriate transformation of the frame. In this 
particular case, it chooses a boost which will cause the |$oo| = | < £ , 22|- 

> T: =dytrgen (Boost , RicciSp) : 

> evalm(T); 

( Tv, v f "i" Tv fv \ - 

W 2 (-rufu+r u , u f)J 



(Tv ,vf~^~ fv \ 
~ .f 2 {-ru.fu + r u , u f)) . 

The dyad transformation T is applied to the curvature spinors using the 
applydytrO function. The following command creates a new frame called 
stdB by applying the spinor dyad transformation, dytr=T, to the curvature 
spinors. (In fact, since WeylSp and Lambda are invariant under T, their 
components are simply copied across to the new frame.) By specifying ro- 
tateTetrad, we also ensure that the frame vectors {I, n, to, to} are themselves 
rotated. 

> applydytr (WeylSp , RicciSp , Lambda , dytr=T , 
rotateTetrad=auto,newname=stdBl) ; 

Transformation completed. 

The rotated frame has been named ' stdBl ' . 

We now switch to the new frame and examine the spinor components. 

> grmetric(stdBl) ; 
Default metric is now stdBl 

> gralter (RicciSp, radsimp) ; 

Component simplification of a GRTensorll object: 
Applying routine radsimp to object RicciSp 

> grdisplay (WeylSp, RicciSp, Lambda) ; 

For the stdBl spacetime : 
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*2 = 

1 94? 1 9e 2 , 1 1 fu, v 1 fu fv 1 9e, e 1 



6 r 2 5 4 6 r 2 g 4 3 r / 3 r 2 / 6 / 2 6 / 3 6 r 2 5 3 6 r 2 g 3 



i \l^v, v I fv \l T u fu ^u, u f 
*0 = 75 

f 2 r 

$! ! = 

i ( 2 g 4 r u r v f + f g^ 2 + f g 2 - g 4 /„, „ r 2 / + .g 4 /„ /„ r 2 - / 3 g g e , e - / 3 5 50, 

/(r 2 .g 4 / 3 ) 



I \j ^v, v f T\> fv \l T u fu Tu, u f 

f r 

. _ 1 30 2 1 g e 2 1 

^ 777 ~~ o a > 777 ~~ o 7 ~r 77 o ^ r 77 ^ r 777 



12 r 2 g 4 12 r 2 5 4 6 r 2 / 3 r / 12 / 2 12 / 3 12 r 2 g 3 
_ _J_ ff0,4> 
12 r 2 5 3 

At this point the frame has been fixed as far as possible, as the curvature 
spinor components are all invariant under the remaining 1-parameter continuous 
isotropy, spins. We can count the number of independent functions among the 
components and subtract the dimension of the frame freedom to determine that 
four scalar polynomial invariants are expected. 

8 Appendix B: Invariants of Class B warped 
product spacetimes 

We record here the invariants {R, r±, r2, w^}- 
For type B\\ 



ds 2 = -2f(«, v) du d v + r(u, v) 2 g(0, (f>) 2 d6 2 + r(u,v) 2 g(6, 



\2 



R = 2(-r(u, v) 2 g{6, 4>f f u f v + r(u, v) 2 g{6, </>) 4 fu, v f(«, v) 

+ 4 r(u, v) g{6, (f>) 4 r u , „ f(«, w) 2 + f(«, v) 3 g^ 2 - f(u, v) 3 ^, g(0, 0) 

+ f(«, V ) 3 3e 2 - f(u, vf g e , e g(0, <P) + 2 f(«, v) 2 r u g(6, </>) 4 r„) / (f(«, «) 3 

r(«, «) 2 g(fl, 0) 4 ) 
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Rl = \ (2 f(«, vf g^ 2 gf + f(«, v) 6 2 g(0, 0) 2 + f(«, vf g e , e 2 g(0, 4>f 

+ 12 r(«, vf g(6, <j>f f u f v f(u, «) 2 r u r„ + r(u, v) 4 g(0, 0) 8 /„, 2 f(u, «) 2 
+ 2 r(u, «) 2 g(0, </>) 5 /„, „ f(«, z;) 4 - 2 r(u, v) 2 g{6, <f>f /„, „ f(u, v) 4 30 2 
+ 2 r(u, «) 2 g(0, 0) 5 /„, „ f(«, vf 98,8-4: r(«, «) 2 g(0, 0) 8 /„, „ f(«, «) 3 r u r„ 

- 2 r(«, w) 2 g(0, 0) 4 /„, „ f(«, vf 30 2 + r(u, vf g(6, <bf f u 2 f v 2 

- 8 f(«, vf r(«, w) 2 g(0, </>) 8 r„, „ /„ r u - 8 f(«, w) 3 r(u, w) 2 g(0, <pf f v r v r u , u 

- 2 f(«, vf 9( f g e , s g(6, 0) + 4 f( U) vf 94> 2 r u g(6, 0) 4 r v 

- 2 f(«, vf 50 , g(0, 0) ge 2 + 2 f(u, w) 6 g<t> , g(0, 0) 2 fle , g 

- 4 f(«, «) 5 5 0, g(0, 0) 5 r u r v - 2 f(u, «) 6 gf ge, e g(0, 4>) 

+ 4 f(«, vf 90 2 r u g(0, </>) 4 r v - 2 r(u, «) 4 g(0, 0) 8 /„ /„ /„, „ f(«, v) 

+ 2 r(u, «) 2 g(0, 0) 4 /„ /„ f(«, v f 5 2 - 2 r(u, «) 2 g(0, 4>f fu U f(u, «) 3 30, 

+ 2 r(u, «) 2 g(0, 0) 4 /„ /„ f(«, vf 9e 2 - 2 r(«, vf g(0, <pf f u f v f(«, v) 3 30, s 

- 4 f(u, «) 5 go, 9 g(0, <pf r u r v +A f(«, w) 4 r„ 2 g(0, 0) 8 r„ 2 

+ 8 f(u, w) 4 r(u, vf g(0, <f>f r v , v r u , u -2 f(u, w) 6 50 2 30, g(0, 0) 

+ f(u, «) 6 30 4 + f(«, «) 6 3/) /(f(«, «) 6 r(«, «) 4 g(0, 0) 8 ) 

= -|(r(u, «) 2 g(0, 0) 4 /„ /„ - r(u, vf g(0, 0) 4 /„, „ f(«, v) + f( U) V ) 3 ff0 2 

- f(«, w) 3 30, g(0, 0) + f(u, vf g e 2 - f(u, w) 3 3e, e g(5, 0) 

+ 2f(u, w) 2 r tt g(0, 0) 4 r t ,)(-r 1 ,, t , f(u, v) + f v r v ) (r u , u f{u, v) - f u r u ) j ( 
f(u, vfr(u, vfg(9, 0) 4 ) 

W2R = u ) 2 g( ' </ > ) 4 /«,,f(w, v) -r(u, vfg(0, <pff u f v 

36 

- 2 r(«, u) g(0, 0) 4 r u , v f(u, w) 2 + f(«, w) 3 50 2 - f(«, «) 3 30, g(0, 0) 
+ f(«, vf 9e 2 - f(«, «) 3 3e, « g(fl, 0) + 2 f(u, «) 2 r u g(6, cj>f r v f /(f(«, «) 9 
r(«, «) 6 g(0, 0) 12 ) 
For type £2: 

ds 2 =f(w, w) 2 rfw 2 +f(u, vf dv 2 -2i(u,vfg(6, (f>) d6 d 4> 

R = 2(r(u, «) 2 g(0, 0) 3 ,/, 2 - r(u, vf g(9, 0) 3 /„, „ f(«, v) + r(u, «) 2 g(0, 0) 3 /„ 2 

- r(u, «) 2 g(0, 0) 3 /„, „ f(«, v) - 2 r(«, v) g(6, <Pf r u , u f{u, vf 

- 2 r(u, v) g(9, (f>f r v , v f(u, w) 2 - f(u, w) 4 30 30 + f(u, w) 4 30, e g(0, 0) 

- f(«, vf r u 2 g(9, 4>f - f( U) vf r 2 g(9, 0) 3 ) /(f(u, «) 4 r(«, «) 2 g(0, 0) 3 ) 
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Rl = 1 (r(u, «) 4 g(9, </>) 6 A, 4 + r(«, vf g(6, <f>f f u 4 + f( U) «) 8 5e 2 ff0 2 

+ f(«, W ) 8 50, e 2 g(0, cj>f + f(«, «) 4 r u 4 g(0, 0) 6 + f( U) «) 4 r„ 4 g(0, 0) 6 

- 2 r(«, w) 2 g(0, 0) 6 /„, „ f(«, vf v 2 - 2 r(u, vf g{9, <pf f v , v f(u, «) 3 r, 2 

- 2 r(«, «) 4 g(0, 4>f f u 2 /„, u f(u, w) + 2 r(«, «) 2 g(0, </>) 3 /„ 2 f(u, «) 4 g e 50 

- 2 r(«, vf g{6, 4>f f 2 f(u, «) 4 50, 9 + 10 r(u, «) 2 g(0, </>) 6 /„ 2 f(u, vf r 2 
+ 10 r(«, v) 2 g(6, 4>f U 2 f(«, «) 2 r, 2 + r(u, «) 4 g(0, «/)) 6 /„, „ 2 f(«, v) 2 

- 2 r(u, v) 2 g(6, 4>f /„, u f(«, «) 5 50 50 + 2 r(u, v) 2 g(0, 0) 4 /„, u f(u, vf 94>t e 

- 2 r(«, vf g(6, 4>f /„, u f(«, «) 3 r„ 2 - 2 r(«, w) 2 g(0, 0) 6 /„, u f(«, «) 3 r, 2 
+ 2 r(u, v) 2 g(0, 0) 6 r„, „ 2 f(«, «) 4 - 8 r(u, vf g(0, 4>f r„, „ f(u, v) 3 /„ r u 

+ 8 r(«, vf g(0, 4>f r u , „ f(«, w) 3 r„ /„ - 4 r(«, «) 2 g(0, 4>f r u , u f(«, «) 4 r„, „ 

+ 8 r(u, V ) 2 g(0, 0) 6 /„ r u f(«, «) 3 r„, „ + 2 r(«, «) 4 g(0, 0) 6 /, 2 /„ 2 

- 2 r(u, vf g(0, 4>f f 2 /„, u f(«, w) + 2 r(«, «) 2 g(0, 0) 3 f 2 f(«, w) 4 50 50 

- 2 r(«, «) 2 g(0, </>) 4 A 2 f(«, v) 4 50, 9 + 10 r(u, «) 2 g(0, 0) 6 /„ 2 f(«, «) 2 r u 2 

+ 10 r(«, z;) 2 g(0, 0) 6 /, 2 f(u, «) 2 r 2 - 16 f(u, v) 3 r(«, «) 2 g(fl, 0) 6 r„, „ /„ r u 

- 16 f(u, vf r(u, vf g(6, 4>f r u , v f u r v -2 f(u, vf g e g^ g</, t 9 g(6, <p) 
+ 2 f(«, vf 9e g<p T U 2 g(9, 0) 3 + 2 f(«, «) 6 3e 94> v 2 g(6, <j>f 

- 2 f(«, vf g^ e g(6, <f>f v u 2 - 2 f(u, vf 9<t>t e g(0, <f>f v 2 

+ 2 f(u, vf r u 2 g(0, r, 2 + 8 f(«, «) 4 r(«, «) 2 g(«, <f>f r u , v 2 

- 2 r(u, vf g(9, <jyf f 2 f v , v f(«, w) + r(«, «) 4 g(fl, 0) 6 /„, , 2 f(u, vf 

+ 2 r(«, vf g(0, 4>f fv, v f(u, «) 2 /„, „ - 2 r(u, «) 2 g(fl, 0) 3 /„, „ f(«, vf g e 50 
+ 2 r(u, w) 2 g(0, 0) 4 /„, „ f(«, «) 5 50, e + 2 r(u, v) 2 g(0, <j)f r v , v 2 f(«, «) 4 

- 8 r(«, «) 2 g(e, 0) 6 r„ /„ f(u, «) 3 r„, „ - 2 r(«, «) 4 g(0, 0) 6 /„, „ f(«, w) / u 2 ) /( 
f(u, «) 8 r( W , «) 4 g(0, 0) 6 ) 

i?^ = --(4f(u, v)f u r u r v , v + Af 2 r v 2 -4f(u, v)r u , u fur u 

-2r u , u f(u, w) 2 r„, „ + r u , „ 2 f(u, w) 2 + 4 / u 2 r M 2 + r^^ t, 2 f(tt, u) 2 + 4 /„ 2 r^ 2 
+ 4/„ 2 r„ 2 + 4r U;t , 2 f(M, w) 2 +4r„,„f(u, u) /„ - 8f(u, v)r u , v f u r v 
-£r v f v {(u, v)r v , v -8r UiV f(u, v)f v r u )(i(u, vfg e g 4> -{{u, vfg^^giO, <f>) 
+ f(«, «) 2 r M 2 g(0, 0) 3 + f( U) «) 2 r 2 g(6, 4>f + r(«, «) 2 g(0, 0) 3 /„ 2 

- r(«, 1;) 2 g(fl, 0) 3 /„, „ f(u, v) + r(«, «) 2 g(0, 0) 3 / u 2 

-r(«, V ) 2 g(tf, ^) 3 /„, u f( W , w)) /(f(«, «) 10 r(u, vfg(6, 0) 3 ) 
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W2R = - gg (r(u, vf g(9, cj>f f 2 - r(u, v) 2 g{6, 0) 3 /„, „ f(u, «) 
+ r(u, v) 2 g(0, </)) 3 / M 2 - r(u, «) 2 g(0, 0) 3 /„, u f(«, v) 

+ t(u, v) g(6, (f>) 3 r„, v f(u, v) 2 + r(u, u) g(0, 0) 3 r u , „ f(w, w) 2 - f(u, u) 4 50 5^ 
+ f(u, «) 4 , g(0, 0) - f(«, «) 2 r u 2 g(0, cf>) 3 - f(«, «) 2 r v 2 g(9, 0) 3 ) 3 / ( 
f(tt, «) 12 r(u, «) 6 g(0, ,/>) 9 ) 

9 Appendix C: Ricci syzygies 

Syzygies for the Ricci invariants (in terms of ri, r 2 and r 3 ) up to rio are given 
by: 

6r 4 - 5rir 2 = 0, (40) 

24r 5 + 3r? - 8r| - 18rir 3 = 0, (41) 

24r 6 - 7r 2 r 1 2 - Ur 3 r 2 = 0, (42) 

144r 7 - 36r 2 - QAr^l - 36r 2 r 3 + 9r 4 = 0, (43) 

36r 8 - 4r| - 27rir 2 r 3 = 0, (44) 

576r 9 - 160r 2 r 3 - WOnrj - 160r 2 r^ + 9r x = 0, (45) 

and 

1728ri - 396r 2 r 2 - 352r ir % - 792r 2 r 2 r 3 + 99r 4 r 2 = 0. (46) 
These syzygies hold in general, unlike the syzygy (|l7|). 
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